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CLOSED  FORM  CHARACTERISTIC  FUNCTION  FOR  GENERAL  COMPLEX 
SECOND-ORDER  FORM  IN  CORRELATED  COMPLEX  GAUSSIAN  RANDOM  VARIABLES 

INTRODUCTION 

The  evaluation  of  the  performance  of  signal  processors  in 
Gaussian  noise  and/or  a  fading  medium  frequently  centers  on 
finding  the  ensemble  average  of  an  exponential  of  a  complex 
second-order  form;  for  example,  see  [1].  Furthermore,  this 
complex  form  can  involve  complex  vectors  with  arbitrary  means 
and  covariance  matrices.  Here,  we  will  address  and  solve  the 
following  specific  problem. 

Random  vector  Z  is  Nxl  with  an  arbitrary  mean  vector  Z  and 

u  np 

covariance  matrices  Z  Z  and  Z  Z  .  Vectors  Dj  and  D2  are  Nxl 
complex  constant  column  vectors,  while  Cj,  C2,  C3,  C4  are  Nxn 
complex  constant  matrices.  The  general  complex  form  c  is  defined 
as 

C  ■  Dj  Z  +  d2  Z*  +  zH  cx  Z  +  zH  c2  Z*  +  zT  c3  Z  +  zT  c4  Z*  .  (1) 

The  quantity  of  interest  is  the  characteristic  function  fc(a), 
defined  here  as  the  ensemble  average  (over  Gaussian  vector  Z), 

fc(a)  =  exp(ac)  ,  (2) 

where  a  is  an  arbitrary  complex  scalar.  This  characteristic 
function  will  be  found  in  its  most  compact  form,  convenient  for 
rapid  evaluation  for  numerous  values  of  a. 
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In  order  to  solve  this  problem,  it  is  necessary  to  solve  a 
generalized  characteristic  value  problem,  with  non-Hermitian 
complex  matrices,  for  the  complex  right  and  left  eigenvectors. 
Accordingly,  the  next  section  is  devoted  to  reviewing  the 
properties  of  generalized  characteristic  value  problems  for 
complex  matrices  and  to  developing  the  necessary  expansions  of 
relevant  determinants  and  matrices.  These  properties  and 
expansions  are  then  used  in  the  succeeding  section  to  derive  a 
closed  form  result  for  characteristic  function  f c ( ex )  in  (2). 
Some  related  work  appears  in  [ 2 ] . 
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RIGHT  AND  LEFT  EIGENVECTORS  OF  GENERALIZED  CHARACTERISTIC 
VALUE  PROBLEM  WITH  NON-HERMITIAN  COMPLEX  MATRICES 

Let  Mxm  matrices  B  and  C  be  complex  non-Hermitian,  and  let 
matrix  C  be  nonsingular.  The  right  Mxl  (column)  eigenvectors 
(RmJ  and  corresponding  eigenvalues  {A^j,  both  of  which  are 
generally  complex,  of  the  generalized  characteristic  value 
problem  satisfy  the  vector  equations 

B  R  *  A  C  R  for  1  1  m  1  M  .  (3) 

m  m  m 

Since  this  can  be  written  as 

(B  -  A  C)  Rm  ■  0  for  1  1  m  £  M  ,  (4) 

the  complex  eigenvalues  { }  can  be  found  from  the  determinant 

det ( B  -  A  C)  =  0  .  (5) 

Because  matrix  C  is  nonsingular,  then  det(C)  ?  0,  meaning  that 
M 

the  term  A  appears  on  the  left  side  of  (5).  Therefore,  there 
are  M  roots  of  (5),  denoted  by  eigenvalues  fAm},  which  are 
assumed  to  be  distinct. 

From  (4),  it  follows  that  the  M  elements  of  the  m-th  right 
column  eigenvector  R^  are  proportional  to  the  M  cofactors  of 
respective  elements  of  any  row  of  the  matrix  B  -  Am  C  (3;  page 
113,  problem  71].  Alternatively,  eigenvector  Rm  is  proportional 
to  any  column  of  the  adjoint  matrix  of  B  -  Am  C.  Thus,  each  Rm 
is  uniquely  determined  by  (4),  except  for  a  complex  scale  factor. 
The  absolute  scaling  of  each  Rm  is  fixed  here  by  requiring  it  to 
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have  unit  length: 

RH  Rm  -  1  for  1  i  m  i  M  .  (6) 

in  in 

However,  each  Rm  is  still  ambiguous  within  a  phasor  exp(i0m), 
where  each  ©m  is  real  but  arbitrary. 

The  set  of  M  right  eigenvectors,  {RmJ,  are  linearly 
independent  when  the  M  eigenvalues  (Am)  are  all  distinct.  To 
prove  this,  assume  that  the  first  n  eigenvectors  are  linearly 
independent,  but  that  the  n+1  eigenvector  is  linearly  dependent; 
that  is,  assume  that  we  can  write 


n 


'n+1  4 — r*  am  Rm  ' 

m=  1 


(7) 


where  at  least  one  of  the  coefficients  lam}  is  nonzero.  Then, 

‘  m 

pre-multiplying  (7)  by  matrix  B,  we  find 


B  Rn+1  *  i  %  B  Rm  ■ 

m=  1 


(8) 


or,  using  (3), 


n 

A  .  C  R_, ,  *  >  ,  am  A_  C  R_  . 
n+ 1  n+ 1  *- — s*  mm  m 

m=  l 


(9) 


Pre-multiply  (9)  by  matrix  C-1,  since  C  is  nonsingular,  to  get 


n 

A,  i  R. .  *  )  *  a  A  R  • 

n+1  n+1  _ — t*  m  m  m 

m=i 


(10) 


Now,  multiply  (7)  by  An+1,  and  subtract  it  from  (10),  to  obtain 

n 


m«l 


(A  -  A  . )  a  R  *  0  . 
'  m  n+l  m  m 


(11) 
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Because  eigenvectors  Rj,...,R  have  been  assumed  linearly 
independent,  all  n  coefficients  in  (11)  must  be  zero.  But,  since 
at  least  one  am  ¥■  0,  then  at  least  one  of  the  eigenvalues 
Xi,...,Xn  must  equal  Xn+^.  This  contradicts  the  assumption  of 
distinct  eigenvalues.  Therefore,  assumption  (7)  is  invalid; 
rather,  all  eigenvectors  are  linearly  independent.  Then,  the 
right  Mxm  normalized  modal  matrix 

R  s  [Rx  •  •  •  Rm]  (12) 

is  nonsingular.  Despite  (6),  it  should  be  noted  that  R  R  ^  I, 
in  general,  due  to  matrices  B  and  C  being  non-Hermitian. 

The  left  Mx 1  (column)  eigenvectors  |Lm}  and  corresponding 
eigenvalues  I Ym 1 /  both  generally  complex,  of  the  generalized 
characteristic  value  problem  satisfy  the  vector  equations 

L®  B  =  Ym  lJJ  C  for  1  £  m  £  M  ,  (13) 

m  m  m 

or 

(B  -  Ym  C)  =  0  for  1  S  m  S  M  .  (14) 

mm' 

The  left  eigenvalues,  lYml»  are  found  from  the  determinant 
equation,  det(B  -  y  C)  =  0,  which  is  identical  in  form  to  (5). 
Therefore,  the  left  and  right  eigenvalues,  of  the  two  generalized 
characteristic  value  problems  (3)  and  (13)  for  matrices  B  and  C, 
are  the  same;  that  is, 

Ym  =  Xm  for  1  £  m  4  M  .  (15) 
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Therefore,  (14)  can  be  written  as 

(B  -  A  C)  =  0  for  1  i  m  i  M  ,  (16) 

m  m  ' 

meaning  that  the  M  elements  of  left  (column)  eigenvector  Lm  are 
proportional  to  the  M  cofactors  of  respective  elements  of  any  row 

LI 

of  the  matrix  (B  -  A  C)  .  Each  Lm  is  uniquely  determined  by 
(16),  except  for  a  complex  scale  factor.  The  absolute  scaling  of 
each  Lm  is  fixed  here  by  requiring  it  to  have  unit  length: 

l"  L  =  1  for  1  <  m  1  M  .  (17) 

However,  each  Lm  is  rtill  ambiguous  within  a  phasor  exp(i«j>m), 
where  each  <|>m  is  real  but  arbitrary. 

The  set  of  M  left  eigenvectors,  { Lm } ,  are  linearly 
independent  when  the  M  eigenvalues  jAmj  are  all  distinct. 

(The  proof  parallels  that  presented  for  the  right  eigenvectors 
in  (7)  -  (12).)  Then,  the  left  Mxm  normalized  modal  matrix 

L  *  [Lx  .  .  .  Lm]  (18) 

is  nonsingular.  Again,  despite  (17),  we  should  observe  that 
H 

L  L  j*  I  in  general,  due  to  matrices  B  and  C  being  non-Hermitian. 
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PROPERTIES  OF  EIGENVECTOR  MATRICES  R  AND  L 

Define  Mxm  complex  diagonal  eigenvalue  matrix 

A  «  diagjAj^  .  .  .  AM]  .  (19) 

Then,  generalized  characteristic  value  problems  (3)  and  (13)  can 
be  written  in  the  compact  matrix  forms 

B  R  *=  C  R  A  ,  LHB«=ALHC,  (20) 

by  means  of  (12),  (15),  and  (18).  All  M  eigenvalues,  { Xm } ,  are 
nonzero  if  and  only  if  matrix  B  is  nonsingular. 

Define  complex  mxm  matrix 

D  =  LH  C  R  .  (21) 

Now,  pre-multiply  the  left  equation  in  (20)  by  LH,  and  post- 
multiply  the  right  equation  in  (20)  by  R.  There  follows, 
respectively,  upon  use  of  (21), 

LHBR=DA,  LH  B  R  «=  A  D  .  (22) 

Since  the  left-hand  sides  are  identical,  we  must  have  D  A  =  A  D. 

However,  because  eigenvalue  matrix  A  in  (19)  is  diagonal,  with 

distinct  element  values,  it  follows  that  square  matrix  D  *  [dmnJ 

must  also  be  diagonal.  The  proof  of  this  very  important  claim  is 

as  follows:  the  m,n  element  of  product  D  A  is  d _ ,  while  the 

r  n  mn 

m,n  element  of  A  D  is  K  d„.  Because  these  elements  must  be 

m  mn 

equal,  this  leads  to  (A„  -  A„.)  d„,„  ■  0.  But,  if  all  the 
1  '  '  n  m '  mn 

eigenvalues  JAm)  are  distinct,  then  we  must  have  dmn  *  0  for 
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m  ^  n;  that  is,  matrix  D  is  diagonal.  However,  matrix  D  is 
generally  not  a  multiple  of  the  identity  matrix  I. 

Denote  the  m-th  diagonal  element  of  D  by  dm;  that  is, 

D  =  diag[dj  .  .  .  dM]  .  (23) 

Given  the  unit-length  eigenvectors  (Rm)  and  {LmJ  (obtained,  for 
example,  via  the  procedures  mentioned  in  the  sequels  to  (5)  and 
(16),  respectively),  we  can  compute,  according  to  (21)  and  (23), 
the  complex  scalars 

dm  =  l”  C  R„  for  1  S  m  S  M  .  (24) 

m  m  m  '  ' 

Then,  from  (21)  -  (24),  we  have  the  two  orthogonality  rules 

lh  C  r  ■  D  -  diagldj  .  .  .  dM]  ,  (25) 

LH  B  R  -  D  A  =  diagldjXj  .  .  .  dMXM]  .  (26) 

All  diagonal  elements  { dm }  of  matrix  D  in  (23)  -  (25)  are 
nonzero.  This  can  be  seen  by  computing  the  determinant  of  (25): 

h  M 

det ( L  )  det(C)  det(R)  -  det(D)  *  TT  dm  .  (27) 

m*l  m 

But,  since  all  three  matrices  on  the  left  side  are  nonsingular, 
the  left  side  must  be  nonzero,  meaning  that  the  product  on  the 
right  side  must  be  nonzero. 

Furthermore,  by  appropriate  choices  of  the  arbitrary  phasors 
{0ml  and  ($mJ  associated  with  the  right  and  left  eigenvectors 
(Rm)  and  {Lm),  respectively,  it  is  easily  arranged  for  all  jdm) 
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calculated  via  (24)  to  be  real  and  positive.  Then,  the  diagonal 
matrix  in  (25)  contains  only  positive  real  elements.  However, 
the  diagonal  matrix  in  (26)  will  still  generally  be  complex, 
because  the  eigenvalues  are  complex. 
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EXPANSIONS  OF  MATRICES 

Various  matrix  expansions  and  determinants  are  possible  from 
the  above  relations.  We  start  with  (20)  and  find,  for  complex 
scalar  0,  that  we  can  express 

C-0B  =  C-  0CRA  R- 1  -  C  R  (I  >  0  A)  R-1  .  (28) 

There  immediately  follows  the  determinant 

M 

det(C  -  0  B)  =  det(C)  det(I  -  0  A)  =  det(C)  TT  (1-0  A  )  .  (29) 

m=l  m 

Once  the  eigenvalues  {Am)  have  been  evaluated  from  (5),  this 
expression  allows  for  rapid  evaluation  of  the  determinant  of 
C  -  (3  B  for  numerous  values  of  complex  scalar  0.  The 
eigenvectors  are  not  required  for  this  calculation. 

From  (21),  we  have  the  useful  result  for  the  inverse, 

R-1  «  D~ 1  LH  C  .  (30) 


When  we  combine  this  with  (28),  we  obtain  a  diagonal  expansion 
for  the  following  matrix: 

(C  -  0  B ) ~ 1  =  R  (I  -  0  A)-1  R~ 1  C_1  =  R  (I  -  0  A)"1  D_1  LH  * 


(1-0 


1 

'•’•'A. 


(31) 


H 

That  is,  terms  of  the  form  R„  L  ,  with  m  /  n,  do  not  appear. 

'  m  n 

Once  the  right  and  left  eigenvectors  have  been  evaluated,  this 
result  allows  for  rapid  evaluation  of  the  inverse  of  matrix 
C  -  0  B  for  numerous  values  of  complex  scalar  0.  Special  cases 
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of  (31)  follow  immediately : 


r- 1  _  1  „  r  H 

C  "  d  Rm  m  ' 

m=l  m 


r2-  1  R  L  H 

r  d  A  m  m 
m=  1  m  m 


(32) 


The  latter  expansion  is  valid  if  all  the  eigenvalues  {AmJ  are 
nonzero,  that  is,  if  matrix  B  is  nonsingular. 

EXAMPLE 

A  numerical  example  of  the  relationships  above,  with  complex 
matrix  elements,  eigenvalues,  and  eigenvectors  is  presented  below 
for  M  =  2.  As  mentioned  above,  the  phasors  associated  with  the 
eigenvectors  have  been  chosen  to  lead  to  a  purely  positive  real 


diagonal  matrix  D. 

The  matrix 

u 

L  R  has  been 

added 

to 

illustrate 

its  nonorthogonal  behavior. 

30-il0 

ilO 

30  - 

15+i5 

B  = 

i 

C  = 

t 

,-18+i6 

8-i6 

-10 

9-i3. 

1 

1  i 

1 

1 

i 

A  =  diag[l+i  1— i ] 

'  R  =  73 

9 

L  =  76 

2+i 

9 

1+i  1+i 

2+i 

.  2+i 

1  3 

r  Ht 

[»  ¥] 

T  Hn  /2 

4+i 

3+i2 

2zi  ! 

3  1 

,  L  L  = 

,  L  R  = 

3 

4+i 

L  C  R  *  D 


10/2 

3 


1  0 
.0  2 


L  B  R  *  A  D 


10/2 

3 


1+i  0  1 

0  2-i2! 


(33) 
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APPROXIMATE  SOLUTION  OF  INTEGRAL  EQUATION 

The  results  above  are  relevant  and  useful  if  we  want  to  solve 
the  following  integral  equation  for  r(y),  with  nonsymmetric 
functions  b(x,y)  and  c(x,y): 

J  dy  b(x,y )  r(y)  -  X  J  dy  c(x,y)  r(y)  ,  (34) 

We  discretize  (34),  leading  to  approximation 

)  '  b(kA,nA)  r(nA)  =  A  )  '  c(kA,nA)  r(nA)  .  (35) 

n  n 

In  terms  of  matrices,  this  is 

B  R  =  A  C  R  ,  (36) 

where  column  vector  R  =  [.  .  .  r(nA)  .  .  .)T  here.  This  equation 
will  have  solutions  only  for  certain  values  of  scalar  A: 

B  R™  *  c  Rm  /  R,®  ( .  .  .  r  ( nA )  .  .  .]T  .  (37) 

m  m  m  m  1  m'  '  J 

At  this  point,  we  have  the  characteristic  value  problem  posed 
earlier  in  (3). 

The  alternative  integral  equation  for  X(x),  where  we 
integrate  on  x  instead  of  y,  is 

J  dx  A*(x)  b(x,y )  ■  y  J  dx  X*(x)  c(x,y)  .  (38) 

This  can  be  discretized,  leading  to  approximation 

5  '  X*(kA)  b(kA,nA)  =  y  HI  X*(kA)  c(kA,nA)  ,  (39) 

k  k 

or,  in  terms  of  matrices. 
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LH  B  =  y  LH  C  , 


(40 


T 

where  column  vector  L  *  [ .  .  .  X(kA)  .  .  .  ]  here.  There  are 
solutions  to  (40)  only  for  certain  values  of  scalar  y: 


L”  B  =  ym  C  ,  L 
m  'mm  m 


Xm(kA)  .  .  .]T  . 


(41) 


This  problem  is  now  identical  to  that  considered  earlier  in  (13). 
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ALTERNATIVE  FORMS  OF  GENERALIZED  CHARACTERISTIC  VALUE  PROBLEMS 

The  generalized  characteristic  value  problems  were  given  in 
(20)  in  the  compact  matrix  forms 

B  R  =  C  R  A  ,  LHB=ALHC.  (42) 

If  we  pre-multiply  the  first  equation  in  (42)  by  C*1,  we  obtain 
an  alternative  matrix  form 

A  R  =  R  A  ,  A  =  C"1  B  ,  (43) 

with  right  eigen-solutions  R  and  A.  The  relation  of  these 
solutions  of  (43),  to  the  earlier  solutions  R  and  A,  is  taken  up 
in  appendix  A.  The  matrix  A  is  not  necessarily  Hermitian,  even 
if  both  matrices  B  and  C  are  Hermitian. 

Similarly,  if  we  post-multiply  the  second  equation  in  (42)  by 
C-1,  we  obtain  yet  another  form, 

LH  A  =  A  LH  ,  A  m  B  C_1  ,  (44) 

with  left  eigen-solutions  L  and  A.  The  relation  of  these 
solutions  of  (44),  to  the  earlier  solutions  L  and  A,  is  also 
addressed  in  appendix  A.  Matrix  A  is  not  necessarily  Hermitian, 
even  if  both  matrices  B  and  C  are  Hermitian.  For  the  example 
given  earlier  in  (33),  the  matrices  defined  in  (43)  and  (44)  are 

Oil  if10  15+i5 

A  =  ,  A  =  i  .  (45) 

-2  2  3 l-3+i  0 

Matrix  A  happens  to  be  real  but  is  nonsymmetric.  Its  eigenvalues 
are  complex,  and  are  given  by  A  in  (33). 
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SPECIAL  CASE  OF  SYMMETRIC  (NOT  HERMITIAN )  MATRICES  B  AND  C 

If  Mxm  complex  matrices  B  and  C  are  both  symmetric/  that  is, 

T  T 

B  =  B  and  C  =  C,  (which  is  not  the  Hermitian  case),  then  (20), 
with  the  left  equation  transposed,  reads  as 

RT  B  ■  A  RT  C  ,  LH  B  ■  A  LH  C  .  (46) 

There  follows  immediately,  for  this  special  case. 


(Actually,  the  phasor  of  each  Lm  could  be  different  from  that  of 
Rm. )  These  relations  allow  for  simplifications  in  the  numerical 
calculation  of  the  complex  eigenvectors;  also,  because  there  is 
now  only  one  fundamental  set  R  (with  L  =  R* ) ,  the  scaling  of  the 
eigenvectors  can  be  altered  to  suit  our  convenience.  In 
particular,  we  now  scale  eigenvector  matrix  R  so  that  diagonal 
matrix  D  in  (25)  becomes  the  identity  matrix.  Then, 
orthogonality  relations  (25)  and  (26)  reduce  to 

RT  C  R  -  I  ,  RT  B  R  -  A  ,  (48) 

while  expansion  (31)  becomes 


(C 


P  B) 


-1 


M 

-  EZ 

m=l 


1  -  P  X 


m 


T 

R  R„ 
m  m 


(49) 
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EXAMPLES  FOR  M  =  2 

The  first  example  has  symmetric  complex  matrices  B  and  C, 
with  the  scaling  of  normalized  modal  matrix  R  now  taken  according 
to  the  first  relation  in  (48). 


1 

-2+i6 

2 

_  1 

i6 

l-i3 

=  2 

2 

-l-i3 

'  C  "  2 

l-i3 

-1  . 

A  =  diag[l+i  1-i]  ,  L  =  R  = 


1 

2 

1+i 

2 


l 

72 

1+i 

72“ 


3 

2+i 

l+i2 

13 

RH  R  = 

4 

2-i 

272 

3 

ii 

K 

E-* 

06 

4 

i3 

272 

-l+i2 

272 

2 

272 

2 

lhcr=rtcr=i  , 


LHBR=RTBR=A=  diag[ 1+i  l-i] 


(50) 


The  only  orthogonality  relations  satisfied  by  the  eigenvectors 
R  and  L  are  the  latter  two,  with  weighting  matrices  C  and  B, 

U 

respectively.  Notice  that  Rm  Rm  ^  1  any  longer;  in  fact,  these 
values  are  3/4,  3/2  for  m  =  1,  2,  respectively. 
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A  second  example,  where  matrix  C  is  now  real  and  symmetric, 
is  presented  next;  matrix  B  is  complex  and  symmetric. 


1  2-i 

3 

1 

B  = 

,  C  = 

2-i  1+i 

1 

2 

=  .83284  -  i. 11405  ,  X2  =  -.63284  +  il. 11405  , 


L  ■  R  = 


,37548  +  i. 07156 
,38980  -  i. 09394 


.51657  +  i. 05202 
-.67808  +  i. 05400 


RH  R  = 


.30687 

06901  +  i. 14124 


06901  -  i. 14124 
.73226 


T 

R  R  = 


27898  +  i. 01949 
07171  +  i. 06008 


07171  +  i. 06008 
,72102  -  i. 01949 


T  -1 
R  R  ■  C 


,4 

,2 


-.2 

.6 


RH  C  R  = 


1.03913 
i. 28249 


-i. 28249 
1.03913 


RTCR*I,  RT  B  R  *  A  »  diag[Xx  X2] 


(51) 
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CHARACTERISTIC  FUNCTION  OF  MOST  GENERAL  COMPLEX 
FORM  WITH  FIRST-ORDER  AND  SECOND-ORDER  TERMS 

Let  Z  be  a  complex  Nxl  random  vector  with  real  and  imaginary 
parts  X  and  Y;  that  is,  Z  =  X  +  i  Y.  Let  D3  and  D2  be  complex 
Nxl  constant  vectors.  Let  C2,  C^,  C4  be  complex  Nxn  constant 
matrices,  which  need  not  be  Hermit ian  or  symmetric. 

The  most  general  first-order  complex  form  in  vector  Z  is 


fl  -  Dj  Z  +  D^  Z*  -  { Dj  +  D2)T  X  +  i(Dx  -  D2)T  Y  «  2  HT  W  ,  (52) 
where 


H 


1 

2 


D1  +  D2 
i(D!  -  D2) 


(53) 


H  is  a  complex  2Nxi  constant  vector  and  is  completely  arbitrary; 
that  is,  every  complex  element  of  H  can  be  independently 
specified.  W  is  a  real  2Nxl  random  vector. 

The  most  general  second-order  complex  form  in  vector  Z  is 


f2  -  ZH  C1  Z  +  ZH  C2  z*  +  ZT  C3  Z  +  ZT  C4  Z*  «  WT  G  W  ,  (54) 

where 


G 


c1  +  c2  +  c3  +  c4 

-i(Cj  +  ^2  -  C3  -  C4) 


i(C.  —  Cn  +  C,  —  C. ) 


C1  “  C2  "  C3  +  C4 


(55) 


G  is  a  complex  2Nx2N  constant  matrix,  which  need  not  be 
Hermitian.  G  is  completely  arbitrary;  that  is,  every  complex 
element  of  G  can  be  independently  specified.  However,  as  far  as 
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the  second-order  form  in  (54)  is  concerned,  only  the  symmetric 

part  of  matrix  G  is  active.  That  is,  we  can  replace  G  in  (54)  by 
.  T 

its  symmetric  part  G  =  (G  +  G  )/2  (which  is  still  complex ) , 

without  affecting  the  value  of  f2;  also,  see  appendix  B. 

T 

Therefore,  we  have  f~  =  W  G„  W. 

2  s 

The  general  complex  form  of  interest  here  is 

c  =  f2  +  fx  =  WT  G  W  +  2  HT  W  =  WT  Gs  W  +  2  HT  W  .  (56) 

The  pertinent  statistics  of  real  random  vector  W  are 

K  =  Cov(W)  =  (W  -  W)  (W  -  W)T  = 


X  -  X 

f  T  -T 

„T  — T 1 

Kxx 

Kxy 

Y  -  Y 

[X  -  X 

Y  -  Y  j  = 

Kyx 

Kyy 

Here,  mean  E  is  a  real  2Nxl  constant  vector,  while  covariance  K 
is  a  real  2NX2N  symmetric  constant  matrix.  Some  important 
aspects  about  obtaining  covariance  matrix  K  from  the  two  Mxm 
complex  covariance  matrices  of  complex  vector  Z  are  discussed  in 
appendix  C. 

We  assume  that  Z  is  a  complex  Gaussian  random  vector.  Then, 

X  and  Y  are  joint-Gaussian,  and  2Nxl  real  random  vector  W  is 
Gaussian  with  probability  density  function 

p(W)  =  (2n)"N  (det  K)~ **  expj-  |(W  -  E)T  K_1  (W  -  E))  .  (58) 

The  statistical  quantity  of  interest  is,  for  complex  scalar  a, 
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the  characteristic  function  of  complex  random  variable  c  in 
(56);  in  particular,  the  characteristic  function  is  defined  here 
as  the  average  of  the  following  exponential: 

fc(a)  =  exp(ac)  =  exp|a  WT  Gg  W  +  a  2  HT  wj  = 

=  J  dW  p(W)  exp(a  WT  Gs  W  +  a  2  HT  w)  = 

=(2n)“N(det  K)-isJdW  exp(-  W-E )TK-1  ( W-E )  +  a  WT  Gs  W  +  a  2  HT  w) 

-  d~h  exp(t)  ,  (59) 

where  (4;  section  8-3] 

d  =  det ( I  -  2  a  G  K)  =  det(I  -  2  a  K  G  )  , 

w  b 

It  -l  -l  lT-i 

t  =  ^  V  ( K  1  -  2  a  Gg )  V-±EAKAE, 

V  =  K-1  E  +  2  a  H  .  (60) 

.  .  T 

Covariance  K  is  a  symmetric  real  2Nx2N  matrix,  Gs  =  (G  +  G  )/2  is 

a  symmetric  complex  2Nx2N  matrix  which  need  not  be  Hermitian, 

mean  E  is  a  real  2Nxl  vector,  H  and  V  are  complex  2Nxl  vectors, 

and  a  is  a  complex  scalar. 

In  general,  we  must  invert  2 N * ? n  real  symmetric  matrix  K. 
Also,  we  must  invert  2NX2N  symmet ( i  complex  matrix  K-1  -  2  a  Gg, 
which  need  not  be  Hermitian,  and  which  depends  on  argument  a. 

If  X  *  0,  Y  »  0,  Dj  «  0,  D2  *  0,  then  E  -  0,  H  ■  0,  V  ■  0, 
t  ■  0,  and  we  need  only  evaluate  complex  determinant 
d  *  det (I  -  2  a  K  Gg),  which  depends  on  a. 
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SPECIAL  CASE  OF  Y  ■  0 


Then,  complex  forms 

fl  =  (D1  +  D2)T  X  =  2  DT  X  , 

f2  =  XT  (C1  +  c2  +  c3  +  C4)  X  ■  XT  C  X  , 

C  =  XT  C  X  +  2  DT  X  =  XT  Cs  X  +  2  DT  X  .  (61) 

Matrices  C  and  D  are  complex  and  completely  arbitrary;  matrix 
T 

C  =  (C  +  C  )/2  is  the  symmetric  part  of  C. 
s 

Identify  in  the  subsection  above, 

2N  -*  N  ,  W  ->  X  ,  E  ->  X  ,  Gs  -»  Cg  ,  H  -*  D  ,  K  Kxx  ,  (62) 

thereby  getting  the  characteristic  function  of  c  as 


(63) 


(64) 


Covariance  K  is  a  real  Nxn  symmetric  matrix,  C  is  a  symmetric 
XX  s 

complex  NxN  matrix  which  need  not  be  Hermitian,  mean  X  is  a  real 
Nxl  vectoi ,  D  and  Vx  are  complex  Nxl  vectors,  and  a  is  a  complex 
scalar. 
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SIMPLIFICATION  OF  GENERAL  CHARACTERISTIC  FUNCTION  (59)  -  (60) 


In  (46),  replace  integer  M  by  2N,  and  replace  symmetric  2Nx2N 
matrices  B  and  C  by  Gg  and  K-*,  respectively.  There  follows 

G  R  =  K_1  R  A  ,  LH  G  =  A  LH  K-1  ,  L  =  R*  .  (65) 

Given  Gg  and  K,  solve  this  characteristic  value  problem  for  2Nx2N 
complex  diagonal  eigenvalue  matrix  A  and  the  corresponding  2Nx2N 
eigenvector  matrix  R.  Then,  orthogonality  relations  (48)  become 

RT  K-1  R  =  I  ,  RT  G  R  -  A  .  (66) 

s 

Now,  in  (29),  replace  p  by  2a,  in  addition  to  the  above. 

There  follows 


-1  -1  1— T 

det(K  -  2  a  Gc)  -  det(K  L)  1 f  (1  -  2  a  A  )  .  (67) 

m=l  m 


Comparison  with  determinant  d  in  (60)  then  yields 


2N 

d  =  det( I  -  2  a  K  Ge)  =  ]“T  (1  -  2  a  A  )  .  (68) 

s  m=l  m 


This  quantity  can  be  quickly  evaluated  for  many  complex  values  of 
a,  once  the  2N  eigenvalues  in  (65)  have  been  determined. 

Also,  from  (49),  we  obtain  the  diagonal  expansion  for  the 
inverse  matrix. 


2  a  Gs) 


1 

2  a  Am 
m 


m 


(69) 


This  enables  the  quantity  t  in  (60)  to  be  expressed  as 
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t  =  l  5ZZ  1 - \ - T“  VT  R_  v  -  \  ET  K-1  E  .  (70) 

2  f — r  1  -  2  a  X  mm  2  '  ' 

m=i  m 


When  we  use  the  definition  of  vector  V  in  (60),  there  follows 


,  2N  (e  +  2  a  h  )  .  „  , 

-T---2  ax7--?  **  «  •  <”> 

m=i  m 


where  we  have  defined  complex  scalars 

em  *  ET  K-1  R„  ,  h  =  HT  R„  ,  for  1  i  m  S  2N  .  (72) 

m  m  m  m 

For  a  =  0  in  (59),  we  have  f  (0)  =  1;  then,  d  in  (60)  or 

v» 

(68)  is  1,  meaning  that  t  must  be  zero  for  a  =  0.  Therefore, 
the  last  term  in  (71)  must  be  expressible  as 


1 

2 


.  2N 
*  m=l 


'm 


(73) 


this  can  also  be  verified  by  use  of  expansion  (69)  for  K  (upon 
setting  a  *  0  there).  Combining  these  results,  (71)  can  be 
simplified  to 


t 


2N  X  e^  +  2  e  h  +  2  a  h^ 
— •  mm  mm  m 


m=l 


1  -  2  a  X, 


m 


(74) 


This  expression  can  be  numerically  evaluated  very  easily  for 
numerous  values  of  complex  scalar  a,  once  the  scalars  in  (72) 
have  been  determined.  Also,  (74)  is  the  most  compact  form  for  t 
and  indicates  the  fundamental  form  for  characteristic  function 
fc(cx)  in  (59). 
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Combined  with  (68),  the  final  result  for  the  characteristic 
function  of  c  is 


fc(ct)  = 


2N 


(m=l 


(i  -  2  a  V 


->! 


exp 


a 


2N  el  +  2  em  hm  +  2  «  hl 

— ■  mm  mm  m 


m=l 


1  -  2  a 


m 


(75) 


This  is  equal  to  [1;  (B-54)]  when  all  the  matrices  are  real.  In 
summary,  (75)  gives  the  result  for  ensemble  average,  exp ( ac ) ,  of 
complex  form 


m  T  T  T  * 

C  -  W  G.  W  +  2  H  w  -  d:  Z  +  Di  2  + 

S  12 

+  ZH  C1  z  +  zH  C2  z*  +  ZT  C3  z  +  ZT  C4  z*  , 


(76) 


where  Z  is  a  random  complex  vector  of  size  Nxl,  { Xm }  and  { Rm }  are 
the  eigenvalues  and  right  eigenvectors,  respectively,  of  2NX2N 
generalized  characteristic  value  problem  (65),  and  coefficients 
{em}  and  {hmJ  are  given  by  (72). 

If  we  expand  the  logarithm  of  this  characteristic  function  in 
(75)  in  a  power  series  in  a,  the  cumulants  (x^(k)}  of  complex 
form  c  are  found  to  be 


XC(D  - 


2N  , 

•  El 

m=l  v 


A.  +  Xe+  2  e  h 
ro  mm  mm 


)  - 


XC(X)  =  2k_1  (k-1) 


2N 


m=l 


,  k-2 
m 


2 

r  +  k  A. 
m 


m  m 


+  h_ 


gi 


for  k  2  2  . 

(77) 


This  agrees  with  [1;  ( B— 6 5 ) ]  when  all  the  matrices  are  real. 
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In  order  to  use  these  results  in  (75)  and  (77),  we  must  solve 
(65)  for  matrices  A  =  diagfAj  .  .  .  X2N]  and  R  *  [Rj  .  .  .  R2n3/ 
and  then  compute  scalars  |em)  and  fhm)  by  means  of  (72). 

Attempts  to  derive  (75)  directly,  by  converting  complex  form 
(56)  to  a  sum  of  (complex)  squares  of  uncorrelated  Gaussian 
random  variables,  have  been  unsuccessful.  We  have  been  unable  to 
discover  any  (complex)  linear  transformation  to  accomplish  this 
goal  directly.  Nevertheless,  we  have  derived  closed  form  result 
(75),  which  is  the  most  compact  form  for  the  characteristic 
function  of  complex  random  variable  c  in  (56).  Additionally,  we 

can  obtain  the  ensemble  average 
^  - -  -  ..  ^ 

exp( 0c* )  =  exp( p*c )  -  (fc(P*)]*. 
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EXAMPLE 

The  exact  result  in  (75)  was  compared  with  simulation 
results  for  the  example  given  earlier  in  (51),  where  M  *  2. 
Specifically,  we  took  N  =  1,  M  =  2N  =  2,  and 


.43 

.77  + 

i .  69 

E  = 

H  = 

• 

(78) 

.31. 

.93  - 

i.  57 

Then,  the  eigenvalues  and 

eigenvectors 

corresponding 

to  matrices 

Gg  and  K-1  are 

=  .8328409  -  i. 1140487  ,  k2  =  -.6328409  +  il. 1140487  , 

[-.3754800  +  i. 0715611  .5165670  +  i. 0520160 

L  =  R  =  .  (79) 

-.3897995  -  i. 0939362  -.6780828  +  i. 0539997 

Notice  that  the  mean  vector  E  and  the  linear  transformation 
vector  H  are  both  nonzero  for  this  numerical  example.  This 
necessitates  the  calculation  and  use  of  eigenvector  matrix  R. 

For  a  =  -1.1  +  il.3,  (75)  yields  the  exact  result 

fc(a)  =  .241135  -  i. 008663  ,  (80) 

while  a  simulation  result  based  on  1E9  trials  of  the  exponential 
of  (56)  gave  the  estimated  value 

f  (a)  £  .241208  ±  .000037  -  i(. 008671  ±  .000028)  .  (81) 

c 
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The  ±  values  are  the  standard  deviations  estimated  during  the 
simulation. 

For  the  alternative  example  a  =  -2.1  +  i3.3,  the  exact  result 
is 

fc(a)  =  .746822  -  i. 410478  ,  (82) 

while  a  simulation  result  based  on  1E9  trials  of  the  exponential 
of  (56)  gave  the  estimated  value 

fc(a)  =  .74619  ±  .00076  -  i(. 41150  ±  .00076)  .  (83) 

One  point  to  observe  about  these  estimates  are  their  rather  large 
standard  deviations,  unless  a  very  large  number  of  independent 
trials  are  taken.  We  will  return  to  this  point  below,  after  we 
have  discussed  the  region  of  convergence  of  the  multiple  integral 
in  ( 59 ) . 
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CONVERGENCE  OF  MULTIPLE  INTEGRAL  IN  (59) 

In  order  that  the  multiple  integral  in  (59)  converge,  certain 
conditions  on  the  complex  2Nx2N  matrix  K_*  -  2  a  Gg  must  be 
satisfied.  To  determine  these  conditions,  consider  general 
complex  2Nx2N  matrix  J  in  the  following  multiple  integral: 

J  dW  exp(-  WT  J  wj  ,  J  =  Jr  +  i  Ji  .  (84) 


Because  only  the  s_  ietric  part  of  matrix  J  is  active  in  ( 84 ) ,  J 
can  be  taken  symmetric  without  loss  of  generality.  Furthermore, 
only  the  real  part  Jr  afiects  the  convergence  of  (84). 

Suppose  we  solve  for  the  2Nx2N  eigenvalue  matrix  y  and 
eigenvector  matrix  Q  of  the  real  symmetric  matrix  Jr: 


Jr  Q  -  Q  Y  ,  QT  Q  -  I  ,  QT  Jr  Q  =  Y 


(85) 


Then,  with  Q  *  [ Q1  .  .  .  Qjjj],  we  have  matrix  expansion 


l  T  2*L 

Jr  -  Q  Y  Q'1  -  Q  Y  QT  -  EZ  Ym  Q*  Qm  . 

m=i 


(86 


allowing  the  real  part  of  the  exponent  in  ( 84 )  to  be  developed  as 


T 

W  Jr  W 


2N,  T  T  2 

I^m  W  °m  %  w  '  £=>m  vin  ' 

m=l  m=l 


(87) 


where  2Nxl  column  vector  V  =  QT  W  »  [v^  .  .  .  v2n^T*  T^e 
transformed  integral  on  V  will  obviously  converge  only  if  ym  >  0 
for  all  m;  that  is,  all  the  eigenvalues  of  real  symmetric  matrix 
must  be  positive.  This  means  that  matrix  Jr  must  be  positive 
definite  in  order  for  integral  (84)  to  converge. 
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When  we  identify  J  with  K~*  -  2  a  G  in  (59),  then 


,-l 


r  =  K  -  2  Re (a  Gg) 


(88) 


For  the  example  in  (78),  this  matrix  is  equal  to 


Jr  = 


3-2a_ 


l-4ar-2ai 


l-4ar-2a^  2-2ar+2a^ 


,  a  =  ar  +  l  ai 


(89) 


This  matrix  is  positive  definite  only  if  ar  <  1.5  and  if 


5  -  2  ar  +  10  ai  -  12  -  4  a?  -  20  ar  ai  >  0  . 


(90) 


These  conditions  on  complex  variable  a  are  depicted  in  figure  1 
as  the  region  inside  the  hyperbola,  which  is  basically  in  the 
second  quadrant.  Any  values  for  a  outside  this  region  lead  to  a 
divergent  integral  (59);  that  is,  the  random  variable  exp(ac)  in 
(59)  has  an  infinite  average  value  when  complex  parameter  a  is 
outside  the  hyperbola  in  figure  1.  In  the  one-dimensional  case, 
this  is  similar  to  the  example  of  a  zero-mean  unit-variance 
Gaussian  random  variable  x  with  transformed  average  value 

exp(ax2)  *  J  dx  (2n)-ls  exp|-  ^  x2  +  a  x2j  (91) 

being  infinite  if  ar  >  1/2. 

In  general,  matrix  Jr  in  (88)  must  be  positive  definite  in 
order  that  (59)  converge,  that  is,  in  order  for  random  variable 
exp(ac)  to  have  a  finite  average.  That  latter  value  is  then 
given  by  (59)  and  (60).  The  region  dictated  by  the  positive 
definiteness  of  Jr  can  take  on  a  variety  of  shapes.  For  example. 
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if  the  2,2  element,  1+i,  of  G$  in  (78)  is  replaced  by  3-i4,  the 
region  of  convergence  of  multiple  integral  (59)  is  the  interior 
of  a  circle  of  radius  3.25  centered  at  a  =  -1.75  -  i2. 5.  Of 
course,  these  regions  always  include  the  origin  a  =  0. 

If  one  tries  to  estimate  the  ensemble  average  of  exp(ac)  by 
an  arithmetic  average  of  a  large  number  cf  independent  trials, 
the  standard  deviation  is  found  to  be  rather  large,  especially 
near  the  boundary  of  convergence  of  integral  (59).  In  fact,  as  a 
approaches  the  boundary,  the  standard  deviation  tends  to  °°.  This 
makes  the  availability  of  closed  form  result  (75)  a  very  useful 
tool  for  analysis,  without  requiring  excessive  computer 
simulation  time. 


Figure  1.  Region  of  Convergence  of  (59)  for  Example  (78) 
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EXTENSION  TO  JOINT  CHARACTERISTIC  FUNCTION 

The  joint  characteristic  function  of  the  real  and  imaginary 
parts,  cr  and  c^,  of  general  complex  form  c  in  (56)  was  derived 
in  [1;  page  141].  Attempts  to  simplify  that  result  have  not  been 
successful.  For  example,  the  determinant  in  [1;  (G-16)  -  (G-18)] 
takes  the  form  det(I  -  l  A  -  C  B),  where  A  and  B  are  constant  Mxm 
matrices,  and  £  and  C,  are  complex  scalars.  Guided  by  (68),  we 
might  hope  that  this  determinant  would  take  on  the  form 

o  M  ,  > 

d  =  det ( I  -  l  A  -  £  B)  =  77  [l  ~  l  -  vm  t)  .  (92) 

However,  this  could  never  be,  because  there  are  not  enough 
degrees  of  freedom  on  the  right-hand  side  of  (92).  To 
demonstrate  this,  consider  the  following  example  for  M  *  2: 

f3  f4  21 

A  -  ,  B  =  ,  (93) 

ll  2  2  3] 

which  are  real  symmetric  matrices.  Then,  the  left  side  of  (92) 
is  given  by 

d  =  1  -  5£  -  7C,  +  5£2  +  13£C  +  8Z,2  .  (94) 

On  the  other  hand,  the  right  side  of  (92)  is  of  the  form 

(1  -  c£  -  Pt)(l  -  /iS  -  v£)  - 

*  1  -  (a+pK  -  (0+\>)£  +  avl2  +  (av+6p)SC  +  •  (95) 

2  2 

If  we  choose  constants  a,/j,0,\>  to  match  the  ,02,  terms,  the 
numerical  values  for  av+Pp  must  then  be  (35±/85)/2  *  12.89  or 
22.11,  which  do  not  equal  the  desired  value  of  13  in  (94). 
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SUMMARY 

The  characteristic  function  of  a  general  complex  second-order 
form  in  correlated  complex  Nxi  Gaussian  random  variables  has  been 
evaluated  in  the  compact  closed  form  (75).  Its  use  requires  the 
solution  of  2Nx2N  generalized  characteristic  value  problem  (65), 
with  complex  non-Hermitian  matrices.  The  final  form  in  (75)  can 
be  accomplished  by  means  of  a  finite  product  and  a  finite  sum  of 
2N  terms,  along  with  one  square  root  and  exponential.  Rapid 
evaluation  for  numerous  values  of  complex  scalar  a  is  readily 
achieved.  The  additional  scalars  required  are  given  in  (72). 
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APPENDIX  A.  ALTERNATIVE  CHARACTERISTIC  VALUE  PROBLEMS 

An  alternative  form  of  the  eigenvalue  problem  was  derived  in 
(43),  namely 

A  R  =  R  A  ,  or  a  R  =  A  R  for  1  S  m  S  M  ,  (A-l) 

—  —  —  —  —  — m  — m  — m 

where  MxM  matrix 

A  =  C~l  B.  ( A-2 ) 

Corresponding  to  these  M  right  eigenvectors  {R^},  there  are  also 
the  left  eigenvectors  { }  of  matrix  A,  which  solve  the  equations 

LH  A  =  v  LH  for  1  i  m  <;  M  .  (A-3) 

— m  —  -Mn  — m 

The  relationship  of  these  solutions  to  the  earlier  results  in 
(3)  and  (13),  for  (Am),  |Rml,  and  { Lm ) ,  is  the  subject  of  this 
appendix.  However,  the  presentation  will  be  significantly 
condensed  here,  because  it  parallels,  to  some  extent,  the 
material  in  the  main  body  of  this  report. 

From  (A-l),  we  find 

(A  -  A*  I)  5*  -  0  ,  det(A  -  A.  I)  -  0  .  (A-4) 

Substituting  (A-2),  this  develops  into 

det(C-1  B  -  A  I)  =  det(C_1 )  det(B  -  A  C)  =  0  .  (A-5) 

'  — m  — m 

But,  since  matrix  C  is  nonsingular,  reference  to  (5)  yields 

A  =  A  for  1  £  m  i  M  ,  ( A— 6 ) 

— m  m 

under  some  ordering  of  the  eigenvalues.  That  is,  the  two 
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eigenvalue  matrices  are  equal:  A  =  A. 
We  can  now  write  (A-l)  in  the  form 


(A  -  \ 


*m  " 


Bm  = 


-  C> 
m 


R  =  0 
—in 


(A-7) 


for  1  £  m  £  M.  When  we  pre-multiply  (A-7)  by  nonsingular  matrix 
C,  and  compare  the  result  with  (4),  we  see  that  we  must  have 


Bm  =  Rm  um  £or  1  *  m  *  M  /  or  R  »  [Rj  .  .  .  1^]  =  R  U  ,  (A-8) 

where  U  s  diag(u^  ...  uM]  is  an  arbitrary  diagonal  matrix  of 
complex  scalar  elements.  The  scaling  of  (F^)  may  be  settled  by 
taking  them  to  be  unit  length;  then,  (6)  indicates  that  we  would 
have  |um|  *  1  for  all  m.  So,  the  right  eigenvectors  {R^  in 
(A-l)  can  be  taken  equal  to  (Rm)  in  (3),  if  desired. 

As  for  the  left  eigenvectors  (L^J,  develop  (A-3)  as 

lJJ  (A  -  Xn  I)  =  0  ,  det(A  -  ^  I)  -  0  .  (A-9) 

Comparison  with  (A-4)  -  (A-6)  yields  =  for  1  S  m  i  M.  We 
can  now  write  (A-9)  in  the  form 

l"  (A  -  k  I)  -  (C_1  B  -  A  I)  =  c'1  (B  -  AC)  =  0.  (A-10) 

— m  m  '  — m  '  m  '  — m  m 

Comparison  with  (16)  yields  the  connection 

c_1  =  v*  ,  or  L  =  CH  L  v  for  1  S  m  i  M  ,  (A-ll) 
— m  m  m  — m  m  m 

using  the  fact  that  C  is  nonsingular,  where  { vm )  are  arbitrary 
complex  scalars.  In  matrix  form,  (A-ll)  is 

L  =  CH  L  V  ,  V  «  diag(Vj  .  .  .  vMJ  .  (A-12) 
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Scalars  (vmJ  in  (A-ll)  can  be  taken  so  that  left  eigenvectors 
|L  1  have  unit  length;  that  is,  |v  I  =  (L^  C  CH  L  )-lj. 

The  relation  for  L  in  (A-12)  can  be  made  to  look  more  akin  to 
R  =  R  U  in  (A-8),  by  using  (25);  namely, 

L-H  =  C-1  L-H  V_H  =  R  D-1  V-H  ■  R  U  ,  ( A- 1 3 ) 

~  —  1  — H 

where  U  =  D  V  is  a  diagonal  matrix. 


PROPERTIES  OF  EIGENVECTOR  MATRICES  R  AND  L 

The  expressions  in  (A-l)  and  (A-3)  can  be  put  into  the  forms 


A  R  =  R  A  ,  LH  A  =  A  LH  ,  (A-14 ) 

with  the  help  of  (A-6).  Now,  define  complex  Mxm  matrix 

D  =  LH  R  .  (A-15) 

Then,  by  arguments  similar  to  those  presented  in  (20)  -  (22), 

LH  A  R  =  D  A  ,  LH  A  R  =  A  D  ,  DA-AD,  (A-16) 
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EXPANSIONS  OF  MATRICES 

A  number  of  useful  expansions  and  relations  can  now  be 
developed  from  the  above  results.  From  (A-14)  and  (A-15), 


M  X14 

An  =  R  An  R_1  =  R  An  D_1  LH  =  >  !  j5  R  lJ  .  (A-20) 

—  —  —  —  —  —  trt  d^  — m  — m  '  ' 

m=  i  — m 

Similarly,  for  complex  scalar  p, 


I  -  P  An  =  R  (I  -  p  An)  R-1  =  R  (I  -  pAn)  D_1LH 


i  1  -  m!  „ 

- A - 2  '  (  A—  21) 

,  d  — m  — m 

m=  1  — m 


n  n  M  n 

det( I  -  (i  aV  det(I  -  P  An)  *  J~1  ( 1  -  P  , 

m=l 


( A-22 ) 


(I  -  p  An)_1  =  R  (I  -  p  An)_1  R-1  =  R  (I 


P  An)_1  D-1  LH  = 


M 


ms 


-  — - - E 

a  d  1  -  p  a" 

— m  m^ 


„  lh 

— m  -m 


(A-23) 


These  diagonal  expansions  allow  rapid  calculation  for  numerous 
values  of  complex  scalar  p,  once  the  eigenvalues  { Am }  and  right 
and  left  eigenvectors  jJ^)  and  jl^}  have  been  determined. 
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EXAMPLE 

For  the  example  considered  earlier  in  (33),  we  have  the 
following  results: 


A  =  C"1  B  = 


0  1 

-2  2 


1 

1  i 

1  ,  1 

1-i  -1+i 

A  =  diag[l+i  1-i]  , 

£  “  73 

1+i  1+i. 

V 

II 

•Jl 

i  1 

RH  R  = 


1 

2+i 

1  "2-il 

3 

_H  _ 

1  3 

2-i 

3 

1 

II 

Jl 

Jl 

-2+i 

[3  1 

L H  R  =  D  =  I  ,  LH  A  R  »  D  A  =  | 


1+i  0 

0  1-i 


( A-24 ) 


The  matrix  A  happens  to  be  real ;  however,  its  eigenvalues  are 
complex.  Relations  (A-8)  and  (A-ll)  have  been  verified  for  this 
example. 
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PROPERTIES  OF  ALTERNATIVE  PRODUCT  A  =  B  C 

Another  form  of  the  eigenvalue  problem  was  stated  in  (44): 

LH  A  =  A  LH  ,  or  LH  A  =  X  LH  for  1  £  m  £  M  ,  ( A-25 ) 

m  m  m 

where  MxM  matrix 

A  =  B  C-1  .  ( A-26 ) 

Corresponding  to  these  left  eigenvectors  {Lm),  there  are  also  the 
right  eigenvectors  ( Rm }  of  matrix  A,  which  solve  the  equations 

A  R  =  y  R  for  1  <  m  £  M  .  (A-27 ) 

m  '  m  m  '  ' 

This  can  be  developed  as  follows: 

(A  -  9m  I)  *  0  ,  det(A  -  Ym  I)  *  0  ,  ( A-28 ) 

det(B  C*1  -  Ym  I)  =  det(B  -  Ym  C)  det(C-1)  =  0  .  (A-29) 

Therefore,  from  (5),  \  for  all  m.  This  enables  the 

'  '  '  m  m 

development  of  (A-28)  according  to 

(A  -  A  I)  R_  =  (B  C'1  -  A_  I )  R  *  (B  -  AC)  C-1  Rm  **  0.  (A-30) 
'  m  m'  m  m  m  m 

Using  (4)  and  the  nonsingular  character  of  matrix  C,  we  find 

C'1  R_  =  R_  u  ,  R=  C  Ru  ,  R  -  C  R  U  ,  (A-31) 

m  m  m  m  mm 

where  U  *  diagfu^  .  .  .  uM] .  This  can  be  compared  with  result 
R  *  R  U  in  (A-8) . 
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Regarding  the  left  eigenvectors  {L^},  (A-25)  yields 

(A  -  X  I)  ■  0  ,  det( A  -  X  I)  =0  .  ( A-32 ) 

m  m  '  '  m 

Comparison  with  (A-28)  -  (A-29)  yields  Xm  =  Am  for  all  m.  Then, 
(A-32)  yields 

(A  -XI)  =  lJJ  (B  C'1  -  A  I)  -  l”  (B  -  Xm  C)  C"1  *  0.  (A- 3  3 ) 
m'  m  '  m  '  m  m'  m 

When  we  post-multiply  (A-33)  by  nonsingular  matrix  C,  we  obtain 


L**  (B  -  A„  C)  =0  . 
m  '  m 


( A-34 ) 


Defining  V  =  diag[Vj  .  .  .  vM],  reference  to  (16)  reveals  that 

( A-35 ) 


rH  _  tH  * 
*m  ^m  vm  ' 


Lm  —  L  v  ,  L  =  L  V  . 
m  mm 


H 


This  can  be  compared  with  result  L  =  C  L  V  in  (A-12). 
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PROPERTIES  OF  EIGENVECTOR  MATRICES  R  AND  L 


The  expressions  in  (A-27)  and  (A-25)  can  be  written  as 


A  R  =  R  A  ,  LH  A  =  A  LH  . 


( A-36 ) 


Now,  define  complex  MxM  matrix 


5  =  lH  R 


(A- 3 7 ) 


Then,  by  arguments  similar  to  those  presented  in  (20)  -  (22), 

LH  A  R  =  6  A  ,  LH  A  R  =  A  5  ,  5  A  =  A  D  ,  ( A-38 ) 


D  *  diagtdj  .  .  .  3MJ  , 


(A- 3 9 ) 


3_  =  Rm  =  um  v*  L“*  C  R  =  u  v*  dm  for  1  £  m  £  M  .  (A-40) 

mmmmmmmmmm 


Thus,  we  have  the  two  orthogonality  relations 


LH  R  =  5  ,  LH  A  R  =  5  A  , 


( A-4 1 ) 


where  both  5  and  A  are  diagonal.  The  elements  {3mJ  can  be 
arranged  to  be  positive  real  for  all  m. 

In  summary,  the  desired  interrelationships  between  the 
alternative  characteristic  value  problems  are 

R=RU,  L  =  CH  L  V  ,  for  A  *  C-1  B  ,  (A-42) 

while 


R*CRU,  L  =  L  V  ,  for  A  *  B  C  . 


( A-43 ) 


Conversely,  having  solved  for  R,L  or  R,L,  these  equations  can  be 
used  to  solve  for  R, L. 
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EXAMPLE 

For  the  example  considered  earlier  in  (33),  we  have  the 
following  results; 


A  *  B  C 


-1  1 


10  15+i5 

•3+i  0 


-  1 

5-i5 

5-i5 

-  l 

1 

-1  1 

'  Rs760 

l+i3 

-3+i 

'  L  =  76 

2+i 

i 

NJ 

1 

H- 

i  5+i 

1  *5~il 

1  T" 

-H  ? 

1  6 

t  L  L  = 

-5+i 

[6  1 

6  1 

rH  s 


/10 

6 


?H  r  * 


L  R  =  D  =  I  ,  L“  A  R  -  D  A  - 


/10 

6 


1+i  0 

0  1-i 


( A-44 ) 


Relations  (A-31)  and  (A-35)  have  been  verified  for  this  example. 
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ALTERNATIVE  FORM  INVOLVING  SQUARE  ROOT  OF  C 

_L 

If  we  pre-multiply  (20)  by  C  ,  and  rearrange  the  equation, 
we  arrive  at 

C-*5  B  C~h  (c*5  r)  -  (c*5  r)  A  ,  (A-45 ) 

which  is  a  more  symmetric  form.  This  leads  us  to  consider  the 
-k 

matrix  C  B  C  ^  and  its  right  and  left  eigen-solutions : 

A  ■  C"1*  B  C-**  ,  A  R  =  R  A  ,  LH  A  =  y  LH  .  (A-46) 

By  a  process  similar  to  that  employed  above,  we  find 

Y  =  A  =  A  ,  R  =  R  U  ,  L  =  L  V  ,  (A-47) 

where  U  and  V  are  arbitrary  diagonal  matrices.  This  can  be 
compared  with  (A-42)  and  (A-43). 
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APPENDIX  B.  GENERAL  SECOND-ORDER  REAL  AND  COMPLEX  FORMS 


Suppose  Kxl  column  vector  X  and  Kxk  arbitrary  matrix  R  are 
real .  Let 


Then, 


Rs  -  j(R  +  RT)  ,  Rd  -  |(R  -  RT) 


T  T 

RZ  •  Re  ,  R j  =  -  R ,  . 

s  s  d  d 


(B-l) 


( B-2 ) 


Consider  the  real  scalar  quantity 


r  =  XT  R  X  =  XT  {Rs  +  Rd)  X  =  rs  rd  , 


where  scalars 


But, 


r  =  XT  R  X  ,  r ,  =  XT  R ,  X  . 
s  s  d  d 


rd  “  rd  *  XT  Rd  X  ■  -  XT  Rd  X  *  -  rd  . 


( B-3 ) 


(B-4) 


( B-5 ) 


Therefore,  rd  =  0,  meaning  that  r  =  rg,  or 


T  T 

r  =  xx  R  X  =  x1  Rr  X  . 

s 


( B-6 ) 


Thus,  we  need  only  consider  the  symmetric  part  of  matrix  R,  at 
least  as  far  as  quadratic  form  r  is  concerned;  that  is,  only  the 
symmetric  part  of  matrix  R  is  active  in  quadratic  form  (B-3). 

More  generally,  suppose  that  Kxl  column  vector  Z  and  Kxk 
arbitrary  matrix  C  are  complex.  Let 


cs  -  i(C  ♦  CH)  ,  Cd  -  |(C  -  CH)  . 


(B-7) 


Then, 


( B-8 ) 
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Now,  consider  the  complex  scalar  quantity 

c  =  ZH  C  Z  =  ZH  (Cg  +  Cd)  Z  =  cg  +  cd  ,  ( B-9 ) 

where  complex  scalars 

cg  =  ZH  Cs  Z  ,  cd  =  ZH  Cd  Z  .  (B-10) 

But, 

cs  =  cs  =  zH  Cs  Z  =  zH  Cs  Z  =  cs  '  (B-ll) 

meaning  that  c  is  real;  also, 

*  H  H  H  H 

Cd  =  Cd  =  ZCd  Z  =  "  ZCd  Z  =  "  Cd  '  ( B— 1 2 ) 

meaning  that  cd  is  purely  imaginary.  Thus,  we  can  express 
(B-9)  as 

c  •  ZH  C  Z  -  cr  +  i  ci  ,  (B-13 ) 

where  the  two  scalars 

cr  =  ZH  Cs  Z  ,  ci  =  -i  ZH  Cd  Z  ,  ( B-14 ) 

are  both  purely  real  quantities.  (If  matrix  C  is  Hermitian,  then 

C  =  C,  Cs  *  C,  Cd  =  0,  and  c^  =  0,  making  c  a  Hermitian  form.) 

It  should  be  noted  that  c,  above  in  »B-9),  is  not  the  most 
general  complex  second-order  form;  tne  latter  is 

C  *  ZH  Cx  Z  +  ZH  C2  Z*  +  ZT  C3  Z  +  ZT  C4  z*  ,  (B-15 ) 

where  C^/  0-2  /  C4  are  arbitrary  complex  Kxk  matrices. 
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APPENDIX  C.  TWO  COMPLEX  COVARIANCES  OF 
A  STATIONARY  COMPLEX  RANDOM  PROCESS 

Let  z(t)  be  a  zero-mean  stationary  complex  random  process 
with  a  complex  covariance,  defined  in  the  usual  way  as 

z(t)  z*(t  -  r )  =  z(t  +  z*[t  -  2  RZZ(T)  *  RZZ("T>  *  (c_1) 

However,  there  is  insufficient  information  in  complex  function 

R„  ft)  to  determine  the  real  auto  and  cross  covariances  of  the 
zz 

stationary  real  and  imaginary  components  of  z(t).  That  is,  with 
z(t)  =  x(t)  +  i  y(t),  (C-l)  yields  the  two  real  equations 

Rxx(T)  +  Ryy(T)  =  Re  Rzz(T)  (even)  , 

Ryx ( T )  -  RXy(T)  "  Im  RZZ(T)  {°dd)  *  (C"2) 

But,  we  cannot  extract  any  of  the  four  real  covariances  on 

the  left  side  of  (C-2),  even  when  we  account  for  the  fact  that 

RyX ( t )  =  Rxy(-r).  Functions  Rxx(t)  and  Ryy(r)  are  always  even 

functions  of  t,  but  R  (t)  need  not  be  odd  in  t.  (As  an  example, 

xy 

consider  y(t)  =  x(t  -  tQ).) 

In  order  to  allow  for  a  unique  solution  of  (C-2),  we  must 
also  know  the  complementary  complex  covariance  of  process  z(t), 

z(t)  z(t  -  T)  -  z(t  +  |)  z(t  -  |)  =  RZ2(T)  =  R„(-T)  •  (C-3) 

(For  a  real  process  z(t),  R^„(r)  and  R__(r)  are  identical.)  The 

z  z  ™*z  z 

complementary  information  in  (C-3)  yields  the  two  additional  real 
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equations 

RXX(T)  -  Ryy ( i )  •  Re  R2Z(t)  (even)  , 

RyX(r )  +  RXy(T)  =  Im  ^ZZ(T)  (even)  .  (C-4) 

The  solutions  to  simultaneous  equations  ( C— 2 )  and  (C-4)  are  now 
available  as 

rkx«>  -  5  Re[Rzz‘T>  +  Ezz'T>]  ' 

Ryy<'>  *  5  -  — z z ( 1 )  ]  ' 

Ryx(t)  =  i  +  SZZ(T>]  , 

RXy(t)  =  i  I">[r«(t)  -  Rzz(t|]  .  (C-5) 

In  the  special  case  where  complementary  covariance  R__(r)  is 

**“Z  z 

identically  zero,  (C-5)  simplifies  to 

Rxx(T)  =  Ryy(T)  *  ?  Re  Rzz(T)  <even>  ' 

Ryx(T)  =  "Rxy(T)  *  \  Im  RZZ(T)  (odd)  *  (C-6) 

When  two  zero-mean  stationary  complex  random  processes  z(t) 
and  w(t)  are  involved,  the  situation  is  slightly  more 
complicated.  We  now  define  the  two  complex  cross  covariances 

z(t)  w*(t  —  x)  s  Rzw(t)  =  R*z(-t)  f 

Z(t)  W(t  -  T)  S  R2W(t)  =  Hwz ( -T  )  .  ( C—  7 ) 
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Then,  if  we  express  the  complex  random  processes  in  terms  of 
their  real  and  imaginary  parts  according  to 


z(t)  =  x(t)  +  i  y(t)  ,  w(t)  =  u(t)  +  i  v(t)  , 


relation  (C-7)  develops  into 


R„  (t)  +  R  (  t  )  +  i  R  (r)-iR  (r)=R  (t), 

xu'  '  yvv  '  yu'  '  xv'  '  zw'  '  ' 


R  (r)  -  R  (t)  +  i  R  ( t )  +  i  R  (t)  =  R  (t)  , 

XU'  #  vv'  7  VII  V  7  W  V  7  — 7U'  7  7 


yv 


yu 


xv  ’ 


-zw' 


the  solutions  of  which  are 


RXU(T)  =  i  Re[R2w(T)  +  RZW(T,J  , 

Ryv(t)  =  5  Re[Rzw<T>  -  Rz„<T>]  < 

Ryu(t)  =  |  I*KW<T)  H-  Rzu(t)]  , 

R*v<r>  -  1  -  RZw«T»]  • 


(C-8 ) 


( C-9 ) 


(C-10) 


The  four  relations  in  (C-10)  must  be  augmented  with  those  in 
(C-5),  as  well  as  with  analogous  ones  giving  the  auto  and  cross 
covariances  of  u(t)  and  v(t)  in  terms  of  Rww(x)  and  R^lx).  The 
end  result  is  that  enough  information  is  available  from  the  six 
complex  covariance  functions  RZ2(t),  R22(t),  Rww(t)f  R^t), 

R  (t),  R„Jr),  to  evaluate  all  of  the  auto  and  cross  covariances 
between  the  four  zero-mean  stationary  real  processes  x(t),  y(t), 
u(t) ,  and  v(t) . 
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EXAMPLE 

Consider  the  pair  of  zero-mean  stationary  complex  processes 
£(t)  and  <a(t)  filtered  by  arbitrary  transfer  functions  A(f)  and 
B(f),  respectively,  as  shown  in  figure  C-l.  The  corresponding 
outputs  are  z(t)  and  w(t),  respectively. 


Figure  C-l.  Two-Channel  Filtering  Operation 

The  two  input  complex  cross  covariances  are 

RCu(T)  *  *-(t)  “*(t  '  T)  '  £c,w(T)  =  «(t  -  t)  .  (C-ll) 

The  two  corresponding  input  complex  cross  spectra  are 

G^(f)  -  J  dt  exp( -i2nfr )  R^w(r)  , 

G^w(f)  -  J  dr  exp(-i2nfr)  F^w(t)  .  ( C- 12 ) 

It  may  be  shown  that  the  two  corresponding  output  cross 
spectra  for  the  complex  processes  in  figure  C-l  satisfy  the  rules 
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Gzw(f)  =  GC,w(f)  A(f)  B*(f)  ' 

Gzw(f)  =  G^w(f)  A(f)  B( -f )  .  (C-13) 

Now,  if  the  product  of  filter  transfer  functions  A(f)  and 
B ( — f )  is  zero  for  all  f,  then  the  complementary  output  cross 
spectrum  G_(f)  is  zero  for  all  f.  This  leads  to  complementary 

™ZW 

output  cross  covariance  fi2W(T)  being  zero  for  all  r,  as  utilized 
in  (C-6).  A  case  where  this  situation  arises  is  when  both 
transfer  functions  A(f)  and  B{f)  correspond  to  one-sided  filters, 
nonzero  only  for  f  >  0.  Processes  z(t)  and  w(t)  in  figure  C-l 
are  then  analytic  processes,  because  auto  spectra  G  (f)  and 
Gww(f>  and  cross  spectrum  Gzw(f)  are  zero  for  f  <  0;  that  is, 
auto  covariances  R  „(t)  and  R.ft)  and  cross  covariance  R„,,(t) 
are  analytic  in  the  upper-half  t -plane. 


51/52 

Reverse  Blank 


TR  10403 


REFERENCES 

[1]  A.  H.  Nuttall,  Exact  Detection  Performance  of  Multiple- 
Pulse  Frequency-Shift  Signals  in  a  Partially-Correlated  Fading 
Medium  with  Generalized  Noncentral  Chi-Squared  Statistics, 
NUWC-NL  Technical  Report  10041,  Naval  Undersea  Warfare  Center, 
New  London,  CT,  23  April  1992. 

[2]  G.  L.  Turin,  "The  Characteristic  Function  of  Hermitian 
Quadratic  Forms  in  Complex  Normal  Variates",  Biometrika,  volume 
47,  pages  199  -  201,  1960. 

[3]  F.  B.  Hildebrand,  Methods  of  Applied  Mathematics, 
Prentice-Hall,  Inc.,  New  York,  NY,  1954. 

[4]  W.  B.  Davenport  and  W.  L.  Root,  Random  Signals  and  Noise, 
McGraw-Hill  Book  Company,  Inc.,  New  York,  NY,  1958. 


53/54 

Reverse  Blank 


TR10403 


INITIAL  DISTRIBUTION  LIST 


Addressee 


Number  of  Copies 


Center  for  Naval  Analyses,  VA  1 

Coast  Guard  Academy,  CT 

J.  Wolcin  1 

Defense  Nuclear  Agency,  VA 

J.  Meyers  1 

Defense  Technical  Information  Center,  VA  12 

National  Radio  Astronomy  Observatory,  VA 

F.  Schwab  1 

National  Security  Agency,  MD 

J.  Maar  1 

National  Technical  Information  Service,  VA  10 

Naval  Air  Systems  Command,  DC  1 

Naval  Air  Warfare  Center,  PA 

Commander  1 

L.  Allen  1 

T.  Madera  1 

A.  Witt  1 

Naval  Command  Control  and  Ocean  Surveillance  Center,  HI  1 

Naval  Command  Control  and  Ocean  Surveillance  Center,  CA 

Commanding  Officer  1 

J.  Alsup  1 

D .  Hanna  1 

W.  Marsh  1 

P.  Nachtigall  1 

C.  Persons  1 

J.  Silva  1 

C .  Tran  1 

Naval  Environmental  Prediction  Research  Facility,  CA  1 

Naval  Intelligence  Command,  DC  1 

Naval  Oceanographic  and  Atmospheric  Research  Laboratory,  CA 

M.  Pastore  1 

Naval  Oceanographic  and  Atmospheric  Research  Laboratory,  MS 

Commanding  Officer  1 

B .  Adams  1 

R.  Fiddler  1 

E.  Franchi  1 

R.  Wagstaff  1 

Naval  Oceanographic  Office,  MS  1 

Naval  Personnel  Research  and  Development  Center,  CA  1 

Naval  Postgraduate  School,  CA 

Superintendent  1 

C.  Therrien  1 

Naval  Research  Laboratory,  FL  1 


D-l 


TR10403 


INITIAL  DISTRIBUTION  LIST  (CONT'D) 


Addressee 


Number  of  Copies 


Naval  Research  Laboratory,  DC 
Commanding  Officer 
D.  Bradley 
W.  Gabriel 
A.  Gerlach 

D.  Steiger 

E.  Wald 
N.  Yen 

Naval  Sea  Systems  Command,  DC 

SEA-00;  -63;  -63D;  -63X;  -92R;  PMS-402 
Naval  Surface  Warfare  Center,  FL 
Commanding  Officer 
E.  Linsenmeyer 
D.  Skinner 

Naval  Surface  Warfare  Center,  MD 
W.  Phillips 

Naval  Surface  Warfare  Center,  MD 
P.  Prendergast 

Naval  Surface  Warfare  Center,  VA 
J.  Gray 

Naval  Surface  Weapons  Center,  FL 
Naval  Surface  Weapons  Center,  MD 
Officer  in  Charge 

M.  Strippling 

Naval  Surface  Weapons  Center,  VA 
Commander 
H.  Crisp 
D.  Phillips 
T.  Ryczek 

Naval  Technical  Intelligence  Center,  DC 
Commanding  Officer 
D.  Rothenberger 

Naval  Undersea  Warfare  Center,  FL 
Officer  in  Charge 
R .  Kennedy 

Naval  Weapons  Center,  CA 

Office  of  the  Chief  of  Naval  Research,  VA 

OCNR-OO ;  -10;  -11;  -12;  -13;  -20;  -21;  -22;  -23  (3) 
P .  Abraham 

N.  Gerr 

D.  Johnson 
A .  Wood 

Space  and  Naval  Warfare  System  Command,  DC 
SPAWAR-00 ;  -04;  -005;  PD-80;  PMW-181 
R.  Cockerill 
L.  Parrish 


1 

1 

1 

1 

1 

1 

1 

6 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

11 

1 

1 

1 

1 

5 

1 

1 


D-2 


TR10403 


INITIAL  DISTRIBUTION  LIST  (CONT'D) 

Addressee  Number  of  Copies 


U.  S.  Air  Force,  Maxwell  Air  Force  Base,  AL 

Air  University  Library  1 

U.  S.  Coast  Guard  Research  and  Development  Center,  CT 

Library  1 

U.  S.  Department  of  Commerce,  CO 

A.  Spaulding  1 

Wright  Patterson  Air  Force  Base,  OH 

R.  Leonard  1 

Brown  University,  RI 

Documents  Library  1 

Drexel  University,  PA 

S.  Kesler  1 

Harvard  University,  MA 

Gordon  McKay  Library  1 

Johns  Hopkins  University,  MD  1 

Kansas  State  University,  KS 

B.  Harms  1 

Lawrence  Livermore  National  Laboratory,  CA 

Director  1 

L.  Ng  1 

Los  Alamos  National  Laboratory,  NM  1 

Marine  Biological  Laboratory,  MA  1 

Massachusetts  Institute  of  Technology,  MA 

Barker  Engineering  Library  1 

Penn  State  University,  PA 

Director  1 

R.  Hettche  1 

E.  Liszka  1 

F .  Symons  1 

Rutgers  University,  NJ 

S.  Orfanidis  1 

San  Diego  State  University,  CA 

F.  Harris  1 

Sandia  National  Laboratory,  NM 

Director  1 

J.  Claasen  1 

Scripps  Institution  of  Oceonography ,  CA  1 

Southeastern  Massachusetts  University,  MA 

C .  Chen  1 

State  University  of  New  York,  NY 

M.  Barkat  1 


D-3 


TR10403 


INITIAL  DISTRIBUTION  LIST  (CONT'D) 


Addressee 


Number  of  Copies 


Syracuse  University,  NY 
D.  Weiner 

United  Engineering  Center,  NY 

Engineering  Societies  Library 
University  of  Connecticut,  CT 
Wilbur  Cross  Library 

C .  Knapp 

University  of  Florida,  FL 

D.  Childers 

University  of  Illinois,  IL 

D.  Jones 

University  of  Michigan,  MI 

Communications  and  Signal  Processing  Laboratory 

W.  Williams 

University  of  Minnesota,  MN 
M.  Kaveh 

University  of  Rhode  Island,  RI 
Library 

G.  Boudreaux-Bartels 
S.  Kay 

University  of  Rochester,  NY 

E.  Titlebaum 

University  of  Southern  California,  CA 
W.  Lindsey 
A.  Polydoros 
University  of  Texas,  TX 
University  of  Washington,  WA 
Applied  Physics  Laboratory 
D.  Lytle 
J.  Ritcey 
R.  Spindel 

Villanova  University,  PA 
M.  Amin 

Woods  Hole  Oceanographic  Institution,  MA 
Yale  University,  CT 

Kline  Science  Library 
A.  Nehorai  * 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


D-4 


TR10403 


INITIAL  DISTRIBUTION  LIST  (CONT'D) 


Addressee 


Number  of  Copies 


Accurate  Automation,  TN 

P.  Simpson  1 

Analysis  and  Technology,  CT 

Library  1 

Analysis  and  Technology,  VA 

D.  Clark  1 

Bell  Communications  Research,  NJ 

J.  Kaiser  1 

D.  Sunday  1 

Berkeley  Research,  CA 

S.  McDonald  1 

Bolt,  Beranek,  and  Newman,  CT 

P.  Cable  1 

Bolt,  Beranek,  and  Newman,  MA 

H.  Gish  1 

EDO  Corporation,  NY 

M.  Blanchard  1 

E  G  &  G,  VA 

D.  Frohman  1 

General  Electric,  MA 

R .  Race  1 

General  Electric,  NJ 

M.  Allen  1 

H .  Urkowitz  1 

Harris  Scientific  Services,  NY 

B.  Harris  1 

Hughes  Aircraft,  CA 

T.  Posch  1 

International  Business  Machines,  VA 

G .  Demuth  1 

Kildare  Corporation,  CT 

R.  Mellen  1 

Lincom  Corporation,  MA 

T.  Schonhoff  1 

Magnavox  Elect.  Systems,  IN 

R.  Kenefic  1 

MSB  Systems,  CT 

A.  Winder  1 

Nichols  Research,  MA 

T.  Marzetta  1 

Northrop  Corporation,  CA 

R.  Nielsen  1 

Orincon  Corporation,  CA 

J.  Young  1 


D-5 


TR10403 


INITIAL  DISTRIBUTION  LIST  (CONT'D) 

Addressee  Number  of  Copies 

1 
1 

1 
1 
1 

1 

1 

1 
1 

1 

1 
1 
1 

1 
1 

1 

1 

1 


D-6 


Prometheus,  RI 
M.  Barrett 

J.  Byrnes 
Raytheon,  RI 

R.  Conner 

S.  Reese 

P.  Baggenstoss 
Rockwell  International,  CA 
D.  Elliott 

Science  Applications  International  Corporation,  CT 

F.  DiNapoli 

Science  Applications  International  Corporation,  VA 
P.  Mikhalevsky 
Sperry  Corporation,  NY 
Toyon  Research,  CA 
M.  Van  Blaricum 
Tracor,  TX 

T.  Leih 
B.  Jones 

K.  Scarbrough 
TRW,  VA 

R.  Prager 

G.  Maher 

Westinghouse  Electric,  MA 
R .  Kennedy 

Westinghouse  Electric,  MD 

H.  Newman 

Westinghouse  Electric,  MD 
R.  Park 


TR10403 


INITIAL  DISTRIBUTION  LIST  (CONT'D) 

Addressee  Number  of  Copies 

Assard,  G.  1 

Bartram,  J.  1 

Beidas,  B.  1 

Bendat,  J.  1 

Breton,  R.  1 

Hahn,  W.  1 

Maltz,  F.  1 

Nicholson,  D.  1 

O'Brien,  W.  1 

Pohler,  R.  1 

Richter,  W.  1 

Schulkin,  M.  1 

Urick,  R.  1 

Von  Winkle,  W.  1 

Werbner,  A.  1 

Wilson,  J.  1 

Admiralty  Underwater  Weapons  Establishment,  Dorset,  England 

Library  1 

Defence  Research  Establishment  Atlantic,  Nova  Scotia,  Canada 

Library  1 

S.  Stergiopoulos  1 

Defence  Research  Estab.  Pacific,  British  Columbia,  Canada 

D .  Thomson  1 

Defence  Scientific  Establishment,  Auckland,  New  Zealand 

L.  Hall  1 

Dept,  of  Science  and  Industrial  Res.,  Lower  Hutt,  New  Zealand 

M.  Poletti  1 

Materials  Research  Laboratory,  Salisbury,  Australia 

Airborne  Sonar  1 

Sonar  Technology  and  Processing  1 

RAN  Research  Laboratory,  Darlinghurst,  Australia  1 

SACLANT  Undersea  Research  Center,  NY 

Library  1 

SIMRAD  SUBSEA,  Horton,  Norway 

E .  Lunde  1 

Sonar  and  Surveillance  Group,  Darlinghurst,  Australia  1 

Welse  Sistemi  Subacquei,  Genova,  Italy 

H.  Van  Asselt  1 


D-7 


TR10403 


INITIAL  DISTRIBUTION  LIST  (CONT'D) 

Addressee  Number  of  Copies 


Canberra  College  of  Advanced  Education,  Belconnen,  Australia 

P.  Morgan  1 

College  of  Technology,  Dublin,  Ireland 

R.  Lawlor  1 

Concordia  University,  Quebec,  Canada 

J.  Krolik  1 

Dalhousie  University,  Nova  Scotia,  Canada 

B.  Ruddick  1 

Indian  Institute  of  Technology,  Madras,  India 

K .  Prabhu  1 

Queensland  University  of  Technology,  Brisbane,  Australia 

B.  Boashash  1 

Royal  Military  College,  Ontario,  Canada 

Y.  Chan  1 

Royal  Naval  Engineering  College,  Plymouth,  England 

I .  Stace  1 

University  of  Alberta,  Edmonton,  Canada 

K.  Yeung  1 

University  of  Auckland,  Auckland,  New  Zealand 

M.  Johns  1 

University  of  Strathclyde,  Glasgow,  Scotland 

T.  Durrani  1 

University  of  Technology,  Leicestershire,  England 

J.  Griffiths  1 

L.  Lloyd,  London,  England  1 

P.  Papoutsanis,  Athens,  Greece  1 


D-8 


